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Abstract 

Detailed numerical studies on Gaussian wavepacket dynamics in quantum systems with symmet- 
ric double- well (DW) potential have been made with the use of the two approaches: the Schrodinger 
approach (SA) and the Heisenberg approach combined with the time-dependent variational approx- 
imation (TDVA) to the Gaussian wavepacket method (GWM). We have calculated time depen- 
dences of wavefunctions, averages of position and momentum, the auto-correlation function and an 
uncertainty product. Our calculations show that a given initial Gaussian wavepacket in the SA is 
quickly deformed at t > where a wavepacket cannot be expressed by a single Gaussian, and that 
assumptions on averages of higher-order fluctuations in TDVA are not justified in the SA. These 
results cast some doubt on an application of TDVA-GWM to double-well systems. 
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I. INTRODUCTION 



The time-dependent state of quantum systems may be described by the Schrodinger or 
Heisenberg approach [Ij. In the Schrodinger approach (SA), the time-dependent wavefunc- 
tion \l/(a;,t) for the one-dimensional system is expressed by the partial differential equation 



ih = H^(x,t) = 

dt ^ ' ^ 



1 & 



,2 



Uix) 



^{x,t), (1) 



2m dx'^ 

where notations are conventional. On the contrary, in the Heisenberg approach (HA), the 
dynamics of the system is described by equations of motion for operators x and p 

d{x) {p) d{p) 



dt m ' dt 



-{U'{x)), (2) 



where the bracket (•) denotes an expectation value. It is generally difficult to obtain exact 
solutions for the Schrodinger and/or Heisenberg equations which are available only for lim- 
ited cases like a harmonic oscillator (HO) system. Much efforts have been made to obtain 
approximate solutions for general quantum systems [Ij. Although equations of motion given 
by Eq. ^ are closed within (x) and (p) for an HO system, they generally yield a hierarchy 
of equations of motion for higher-order fluctuations such as {Sp'^) and {6x6p + 6p6x) 

where 6x = x — {x) and 6p = p — {p). 

The Gaussian wavepacket method (GWM) is one of such methods whose main aim is a 
semi-classical description of quantum systems (for a recent review on GWM, see Ref. [2j). 
If the wavefuction is Gaussian at t = in an HO system, it remains at all t > 0. Heller [3] 
proposed to adopt a (thawed) Gaussian wavepacket given by 



^^if(x,t) = exp 



^[A{x-{x)f + {p){x-{x))+^] 



(3) 



even for more realistic potentials, where A and 7 are time-dependent complex parameters. 
Heller [3j derived equations of motion for (x), (p), A and 7, employing an assumption that 
the potential expanded in the Taylor series at x = {x) may be truncated by 

U{x) = U^'\{x)) + U^'\{x)){x - {x)) + lu^^\{x)){x - {x)f. (4) 

The concept of the Gaussian wavepacket has been adopted in many flelds [2]. Dynamics is 
well described by the GWM for an HO system where motions of fluctuations are separated 
from those of {p) and (g), and they lead to the uncertainty relation: {5q^){5p^) > h^/A. 



Various types of variants such as the frozen [4] and generahzed Gaussian wavepacket methods 
[5] have been proposed [2j. Among them, we pay our attention to the time-dependent 
variational approximation (TDVA) which employs the normalized squeezed coherent-state 
Gaussian wavepacket given by |£,-JJ)J 



(5) 



4/i ' ' h 

fi and a being time-dependent parameters. For the introduced squeezed coherent state, 
Heisenberg equations of motion given by Eq. ^ are closed within (x), (p), {Sx'^) and 
{5x5p + 5p5x) [see Eqs. (33)-(36)]. A comparison between Heller's GWM and the TDVA is 
made in Refs. P[l0]. 

There have been many studies on the GWP (including TDVA) which is applied to HO, 
anharmonic oscillator (AO) and Morse potentials [2j. The GWM has, however, some diffi- 
culty when applied to potentials including terms of x^ with n > 2. Although it has been 
claimed that the GWM yields a fairly good result for AO systems [6j, we wonder whether 
it actually works for double-well (DW) systems. DW potential models have been employed 
in a wide range of fields including physics, chemistry and biology (for a recent review on 
DW systems, see Ref. [llj). Lin and Ballentine [12j, and Utermann, Dittrich and Hanggi 
[13] studied semi-classical properties of DW systems subjected to periodic external forces, 
calculating the Husimi function [14j. They show a chaotic behavior, in accordance with 
classical driven DW systems. Igarashi and Yamada [15] studied a coherent oscillation and 
decoherence induced by applied polychromatic forces in quantum DW system. By using the 
TDVA, Pattanayak and Schieve [8j pointed out that a chaos is induced by quantum noise in 
DW systems without external forces although classical counterparts are regular. This is in 
contrast to the usual expectation that quantum effects suppress classical chaos. Chaotic-like 
behavior was reported in a square DW system obtained by the exact SA calculation [16j. 
The validity of quantum chaos pointed out in Ref. is still controversial [IT7H22]. 

The purpose of the present paper is twofold: to study dynamics of a given initial Gaussian 
wavepacket in DW systems with the use of both the SA and the HA combined with the 
TDVA, and to make a comparison between results of the two approaches. Such a study has 
not been reported as far as we aware of, and it is important to clarify significance of the 
TDVA-GWM applied to DW systems. 

The paper is organized as follows. In Sec. II, we will mention the calculation method 
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employed in our study. We consider DW systems described by the symmetric DW (SDW) 
model. In solving the time-dependent Schrodinger equation, we have adopted the spectral 
method with energy matrix elements evaluated for a finite size A^^ (= 20 and 30). Our 
method is applied to SDW model in Sec. Ill, where we report time-dependences of the 
magnitude of wavefunction an expectation value of x ((x)), the uncertainty 

product {{8x^){5p^)) and the auto-correlation function {C{t)). In Sec. IV we apply our 
method also to an AO model. Sec. V is devoted to our conclusion. 



II. THE ADOPTED METHOD 
A. Schrodinger approach 

We consider a quantum system whose Hamiltonian is given by 

i/=|^ + C/W, (6) 

where m, x and p express a mass, position and momentum, respectively, of a particle and 
U (x) stands for the DW potential. Various approximate analytical and numerical methods 
have been proposed to solve the Schrodinger equation given by Eq. ([T]) [Ij. In our study we 
solve the Shrodinger equation, treating the Hamiltonian H as 

H = Ho + U{x) - Uo{x) = Ho + V{x), (7) 

with 

Ho = ^ + Uo{x), (8) 

V{x) = U{x)-Uq{x). (9) 

Here Hq denotes the Hamiltonian with an appropriate potential [/q, for which eigenfunction 
and eigenvalue are assumed to be given by (l)n{x) and Eq^^ respectively. 

For the stationary state, we expand the eigenfunction ^^(x) in terms of (/)n{x) 

Nm 

^{X) = ^CnMx), (10) 

leading to the secular equation 

n=0 



where E denotes the eigenvalue and A^^ is the maximum quantum number. 

For the time-dependent state, we adopt a spectral method, expanding the eigenfunction 
in terms of 

Nm 

= Y^Cn{t)(l)n{x). (12) 

Time-dependent expansion coefficients {cn{t)} obey equations of motion given by 

Or (f) 

in = Y^^rnn Cn{t) (m = to Nj, (13) 

n=0 



with 



/oo 
V{x)(l)n{x) dx. (14) 
-OO 



Equation (13) expresses the (A^^ + 1) first-order differential equation, which may be solved 
for a given initial condition of {c^(0)}. Initial values of expansion coefficients {c^(0)} are 
determined by 

/oo 
MxY ^g(x, 0) dx (n = to N^), (15) 
-OO 

for a given initial Gaussian wavepacket 0) 

(16) 



*G(-,0) = ^^^exp 



[l-ia) 2 , .Po(^-^o) 



4/x h 
where ji and a are assumed parameters, and and are initial position and momentum. 



respectively, at t = 0. Once solutions of {cn{t)} in Eq. (13) are obtained, the wavefunction 
^{x^t) may be constructed by Eq. ( [T2| ). 

In an alternative spectral method, the solution of the time-dependent Schrodinger equa- 
tion may be expressed by 



Nrr 



= J^a,^,(x)e-^^^^/^ (17) 

where ^^^{x) and E^, are eigenfunction and eigenvalue of the stationary state given by 

H^,{x) = E,^,{x) (z/ = to NJ. (18) 

The time-independent expansion coefficient is determined by a given initial Gaussian 
wavepacket 

/oo 
"^uixY ^^g(^,0) dx. (19) 
-oo 



Both spectral methods given by Eqs. (12) and (17) yield the same results. The advantage 
of the former spectral method is that we may easily calculated the auto-correlation function 



and expectation values expressed in terms of {cn{t)} [see Eq. (54)- (59)], while the latter 
spectral method has clearer physical meaning. Model calculations to be reported in the 
following Sec. Ill have been made by using the former spectral method. Wavefunctions 
obtained by the spectral method have been cross-checked by the MATHEMATICA resolver 
for the partial differential equation. 



B. Heisenberg approach 



It is possible to study dynamics of the quantum system with the use of the HA. Heisenberg 
equations of motion are expressed by P4T0] 

d{x) _ (p) 



dt 
dt 

d{5x^) 
dt 

d{5x5p + 5p5x) 
Jt 



m 



k=2 



k\ 



m 



-{6xSp + 6p6x)^ 



k=i 

oo 



m 



k=i 



k\ 



(20) 
(21) 

(22) 
(23) 

(24) 



dt 

where Sp = p — {p) and 6x = x 
which are not closed in general. 

It is possible to construct various approximations depending on how many terms are 



x). Equations (20)- (24) include high-order fluctuations 



taken into account in Eqs. (20)- (24). If we neglect the second term of Eq. (21), Eqs 



(20) and (21) form classical equations of motion. When we neglect the second term in Eq. 



(21) and truncate Eqs. (23) and (24) at = 1, Eqs. (20)-(24) reduce to equations of 



motion in Heller's Gaussian wavepacket approximation. Equations of motion including up 
to fourth-order corrections were obtained in Ref. fTO]. 

One of conceivable ways to close equations of motion is to assume that a wavepacket is 
expressed by the normalized squeezed coherent state given by Eq. g [6l|8Hl0], which leads 
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to the relations 



{5x5p + 5p5x) 



II {Sx^'^^)=0, a = 1,2, 



4:fl ' 



a. 



(25) 

(26) 
(27) 



/i and a being time-dependent parameters in the wavepacket. Equations (25)- (27) yield the 
uncertainty product 

h"^ + {5x5p + 5p5xY 



{Sx'){Sp') 



From these relations, equations of motion given by Eqs. (|20|)-(|24|) become 

d{x) 



dt 

d{p) 
dt 



M 

m 



da + a'^ 



dt 



2fi 



^ U(^'^H{x)) 



(^-1)! 2^- 



Alternatively, Eqs. (29)- (32) are rewritten as 



d{x) 

~dr 

d{p) 
dt 

d{Sx^) 
dt 

d{5x5p + 5p5x) 
Jt 



m 



2\i 



i=l 



e\ 2^ 



1 



m 



{5x5p + 5p5x)^ 



+ {5x5p + 5p5x) 
2m{6x'^) 



f/(2^)((x))(5x2)^ 



^ (£-l)!2^- 



which are closed for (x), (p), {5x'^) and {5x5p^5p5x). 
By using a change of variables given by [6l IH [9j 



2p 



(28) 



(29) 
(30) 

(31) 
(32) 



(33) 
(34) 

(35) 
(36) 



(37) 
(38) 
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Eqs. (29)-(32) are transformed to 

d{q) (p) 



dt 

dt 

dp 
dt 



m 



dt Ap^ (£- 1)! 2^-1^^ ■ 



(39) 
(40) 

(41) 
(42) 



It was shown that fluctuation variables p and tt are conjugate and that the effective Hamil- 
tonian may be expressed in the extended phase space spanned by {p) , ((?) , p and tt [6l [H |9] 



eff 



2 TT^ 



2m 2m Smp^ 



-2e 



(43) 



We should note that the effective Hamiltonian given by Eq. (43) rehes on identities given 



by Eqs. (25)- (27) which are based on the assumed squeezed Gaussian wavepacket given by 
Eq. ([5]). If they are not held, the effective Hamiltonian given by Eq. (43) is not valid in 
DW systems: related discussion is given in Sec. V. 



III. MODEL CALCULATIONS 

We apply our calculation method to the symmetric DW potential given by 

Uix) = (^) {x'-xir, (44) 

which has stable minima at x = ±Xs and an unstable maximum at x^^ = with the potential 
barrier of A = U{0) — U{±Xs) = muj'^xl/S. As for Uo{x)^ we adopt an HO potential given 

by 

Uo{^) = , (45) 

where uj denotes the oscillator frequency. Eigenfunction and eigenvalue for Hq are given by 



r 



Eon = {n+l]hw (n = 0,l,2-,---), (47) 



where Hn{x) stands for the Hermite polynomials. 



Matrix elements Hjnn in Eq. (14) are given as follows: We rewrite the potential U{x) as 

A^x"^ Asx^ A2X^ 



U{x) 



4 



+ Aix + Ao, 



with 



A, 



mu 
~2¥ 



, Ai = 0, Ao = 



, ^3 = 0, A2 



By using relations given by 



^(at + a), p = i — {a^-a), {9 = ^) 



a) (pn = Vn+1 (pn+l, Cl4>n = Vn4>n-1, 



(48) 



(49) 



(50) 
(51) 



we obtain the symmetric matrix elements H^n for m < n given by 



{n + l/2)hw + 



SA^ 
16 



(2n2 + 2n + 1) + ^(n + 1/2) + Aq 



+ 



^m,n— 1 



A', 9 



{n - l)^yn{n - 1) + -^^/n{n - 1) 



^m,n—2 



+ ^ (I) ^ ^"-(77- - l)(n - 2) 5^,n-3 

H ^\/n{n - l)(n - 2)(n - 3) 6m,n-i, 

Id 



(52) 



where ^2 = A2 — muj^ . 

Various time-dependent quantities may be expressed in terms of {cn{t)} as follows: After 
some manipulations with the use of the relations given by Eqs.([50|) and ([51]), the auto- 
correlation function is given by 



C{t) 



'^{x,t)*'^{x,{)) dx, 



N„ 



= 5]c„(t)*c„(0), 



(53) 
(54) 



n=0 
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and expectation values of x{t) and p{t) are expressed by 



{x{t)) 



J2 4^l{t)Cn{t) + V^4_i(t)c„(t) 



'^2 



(pit)) = ^V^I][V^<+i(t)c„(t)-V^<_i(t)c„(t) 



(55) 
(56) 



{<tf) = (I) J2^Vin + l)(n + 2) <+2(t)c„(t) + (2n + 1) <(t)c„(t) 

+ V^(n-l)<_2(t)c„(t)], (57) 
J][V(n + l)(n + 2)<+2Wcn(t) - (2n + 1) <(t)c„(t) 



+ V^(n-l)4_2(t)c„(t)], 
(a;(t)p(a;)+p(t)x(t)) = i + + 2) <+2(t)c„(t) 



(58) 



- V^(n-l)<_2(t)c„(t)]. (59) 
Figure hi expresses the adopted quartic DW potential U{x) with Xg = and A = 1.0 




FIG. 1: (Color online) The symmetric DW potential (solid curve) with Xs — 2\/2 and A = 1.0 in 
Eq. (44), dashed curves expressing eigenvalues ol (z/ = — 4). 
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FIG. 2: (Color online) The classical x-p phase space of the SDW model for various initial energies 
{E(0)), marks x denoting positions of initial states adopted in our calculations. 



in Eq. (44). We have calculated energy matrix elements of i/rnn, by using Eq. (52) with 
m = 1.0, uj = 1.0, h = 1.0 and A^^ = 30. Obtained eigenvalues are = 0.450203, 0.474126, 
1.09262, 1.39334 and 1.91286 for u = to 4, respectively, which are plotted by dashed curves 
in Fig. [l} The ground state (£"0) and first excited state (£"1), which are quasi-degenerate, are 
below the potential barrier of A = 1.0. The energy gap between £"0 and Ei is 6 = 0.023923. 
Low-lying eigenvalues calculated with A^^ = 30 are in good agreement with those obtained 
with A^^ = 20. 

Figure [2] shows the classical x-p phase space for initial energies of E{0) = 0.0, 1.0 and 
1.125. Marks x in Fig. [2] show two kinds of initial states in the x — p phase space: (1) 
(xo,Po) = (—2^/2,0.0) and (2) (xo,Po) = (0.0,0.5). Calculated results for the two initial 
states (1) and (2) will be separately reported in the following. 
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1. Case of the initial state of (xo,po) = (— 2\/2, 0.0) 



We have adopted the initial Gaussian wavepacket ^g(^) locating at the stable point of 
the left well with (xo,Po) = (—2^/2,0.0), ji = 0.1 and a = 0.0, which yields the mini- 
mum uncertainty product of {5x^){5p^) = 1/4 at t = 0.0. A norm of the initial Gaussian 



wavepacket is c^(0)*c^(0) = 0.999999. Initial coefficients {c^(0)} calculated by Eq. (15) 
are real with appreciable magnitudes at 3 < n < 10. After solving (A^^ + 1) first-order 



differential equations for {cn{t)} given by Eq. (13) for initial values of {c^(0)}, we obtain 
the time-dependent eigenfunction expressed in terms of {cn{t)} in Eq. ( [T2| ). 

Figure [3] shows the 3D plot of calculated by the SA. We note that the initial 

Gaussian wavepacket in the SA is quickly spread as the time develops. In order to scrutinize 
the behavior of t) p at small t, its time dependence at < t < 25 is plotted by bold solid 
curves in Fig. [4| where solid curves denote results of the TDVA. The Gaussian wavepacket 
at t = 0.0 in the SA is widespread even at t = 5.0, and its trend becomes more significant 
with increasing t. Figure [3] clearly shows that in the SA is quite different from 

that in the TDVA and that cannot be expressed by a single Gaussian wavepacket 

except at t = 0.0. 

The difference between the SA and TDVA is more clearly seen in the time-dependent 
expectation value of {x). Figure [5](a) shows {x) of the SA, which expresses a tunneling 
of a particle oscillating with the period of about 260, which is consistent with the period 
estimated from the energy gap by T = 27t/6 = 262. On the contrary, {x) of the TDVA in 
Figs. [5](b) shows more rapid oscillation with a period of about 25 — 30. 

Figures |6](a) and|6](b) show (x) vs. {p) plots calculated by the SA and TDVA, respectively. 
The (x) vs. (p) plot of the SA in Fig. [6](a) is quite different from that of the TDVA in Fig. 
lb). 

Figure [7|(a) shows the uncertain product of {5x'^){5p'^) calculated by the SA, which ex- 
presses a measure of quantum fiuctuation. It starts from the minimum uncertainty of 1/4 at 
t = 0, and with increasing t it grows and oscillates between about 5 and 17 with the period 
of about 130. For a comparison, we plot {l + {5x5p + 5x5py)/4 in Fig. [7|(b). In the TDVA, 
both the quantities agree as expressed in Eq. ([28]). Figures [7|( a) and [7|(b) show that this 
equality is not satisfied in the SA. 

Figures [8](a) and |8](a) show the auto-correlation functions calculated by the SA 
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FIG. 3: (Color online) 3D plot of as functions of x and t of the symmetric model 

(xo = -2Apo = 0.0). 



and the TDVA, respectively, with Eq. ([54|). \C{t)\^ of the SA which is unity at t = 0.0, 
oscillates between about 0.1 and 0.7 with a period of about 260, same as that of (x) shown 
in Fig. [5](a). The result of the TDVA is again rather different from that of the SA. 

2. Case of the initial state of (xo,po) — (0.0, 0, 5) 

Next we adopt a different initial Gaussian wavepacket locating at (xo,Po) = (0.0,0.5) 
with /X = 0.1 and a = 0.0, for which the initial energy is £"(0) = 1.125 (see Fig. [2]). Initial 



coefficients {c^(0)} calculated by Eq. (15) are complex with appreciable magnitudes at 
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FIG. 4: (Color online) The x dependence of at various t of the SDW model calculated 

by the SA (bold solid curves) and the TDVA (solid curve) {xq = — 2\/2,po — 0.0). 

^ ^ 15. The initial state of (xo,Po) = (0.0,0,5) locates near a top of the potential 
barrier (see Fig. [2]). Note that in the classical calculation, the x vs. p plot forms a cocoon 
shape extending from x = —4.06021 to x = 4.06021 and from p = —1.5 to 1.5, as shown in 
Fig. [2| Then at t > 0, a classical particle rolls down the potential up to x = 4.06021 and 
approaches x = —4.06021 after passing through x = in the classical calculation. However, 
this classical behavior is quite different from quantum results calculated by the SA and 
TDVA. The 3D plot of |^^(x, t)p of the SA shown in Fig. [9] has an appreciable magnitude at 
~5 ^ X < 5. Bold solid curves and solid curves in Fig. [To|show |^^(x,t)p calculated by the 
SA and the TDVA, respectively. |^^(x,t)p in the SA is distorted and widespreads at t > 0, 
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FIG. 5: (Color online) The time dependence of (x) of the symmetric model calculated by (a) the 
SA and (b) the TDVA (xq = -2\/2,po = 0.0). 

which is different from the relevant result of the TDVA. An expectation value of (x) of the 
SA in Fig. [TT](a) does not so much depart from the initial point of x = 0.0 in contrast to 
that of the TDVA shown in Fig. [lT|(b). 

The (x) vs. {p) plots calculated by the SA and TDVA are shown in Figs. [l2](a) and 
[T2|(a), respectively. The result of the SA in Fig. [T2|(a) shows a chaotic-like motion, which is 
different from a quasi-periodic motion of the TDVA in Fig. [T2|(b). 



Figure [l3|(a) and 13[h) shows SA calculations of {8x^){5p^) and (1 + {8x5p + 5x5]9)^)/4, 
respectively. The TDVA yields that the former quantity is equal to the latter, which is not 
realized in the SA. 
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FIG. 6: (Color online) The (x) vs. (p) plots of the symmetric model calculated by (a) the SA and 
(b) the TDVA, time step being At = 1.0 (xq = -2\/2,po = 0.0). 

IV. DISCUSSION 

We have studied Gaussian wavepacket dynamics of quantum DW systems in the preceding 
section. It is worthwhile to examine also an AO model given by 

U{x) = y + ^ = f/o(x) + ^, (60) 

where h expresses a degree of anharmonicity. We have repeated calculations, by using the 
method mentioned in Sec. II with necessary modifications. 

Figure [T4|(a) shows the 3D plot of \^{x^ t)\^ for h = 0.01 with an assumed initial Gaussian 
state for (xo,Po) = (—1.0,0.0), ji = 0.1 and a = 0.0. In the case of the HO {b = 0.0), 
|^^(x,t)p is periodic with a period of Tq = 27r/uj. Figure [l4|(a) shows that |^^(x,t)p for 
a small b = 0.01 is nearly periodic as that for b = 0.0 at t < 10. For a larger 6 = 0.1, 
however, this periodicity is destroyed and the wavepacket spreads in the non-Gaussian form 



as Fig. 14:[h) shows. This is more clearly realized in the x-dependence of |^^(x,t)p for the 



SA plotted by bold solid curves in Fig. [T5| which are diflFerent from those of the TDVA 
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FIG. 7: (Color online) Time dependence of (a) the uncertainty product: {5x^){5p^) and (b) 
(1 + {5x5p + 5x5p)'^)/A of the symmetric model calculated by the SA (xq = — 2\/2,po = 0.0). Note 



that {6x^){6p^) equals to (1 + {6x6p + 6x6p)^)/4: in the TDVA [Eq. (28)], which is not realized in 
(a) and (b). 



shown by solid curves. Our result of the SA in Fig. 15 is consistent with that in Ref. 
which studies effects of anharmonicity and interactions in DW systems. 

Figure [l6](a) and [l6](a) show time dependences of (x) calculated by the SA and TDVA, 
respectively, {x) oscillates with a period of about 5.6 in both results. However, a period of 
its envelope variation in the SA 110) is larger than that in the TDVA 30): the former 
corresponds to the revival time after which a wavepacket periodically returns to the initial 
shape. 

V. CONCLUSION 

By using the SA and HA combined with the TDVA, we have numerically studied dynamics 
of Gaussian wavepackets in quantum DW systems with the symmetric quart ic potential. A 



17 




FIG. 8: (Color online) Time dependences of the auto-correlation function of the symmetric 

DW model calculated by (a) the SA and (b) the TDVA {xq = -2\/2,po = 0.0). 

comparison between calculated results of the SA and TDVA-GWM has shown the following: 

(1) With time developing, a given initial Gaussian wavepacket of the SA spreads and deforms 
in symmetric DW system where a wavepacket at t > cannot be expressed by a single 
Gaussian, 

(2) Time dependences of expectation values of (x) and (p) and the auto-correlation function 
C{t) in the TDVA are rather different from their counterparts in the SA, and 



(3) The relation for uncertainty product given by Eq. (28) derived in the TDVA is not 
satisfied in the SA. 

The item (1) which holds also in asymmetric DW systems [26j, is in contrast to the postulate 
expressed by Eq. ^ or Eq. ^ in the GWM. Items (2) and (3) suggest that TDVA cannot 
be a good approximation for DW systems, which throws doubt on studies having been 
made with the use of the TDVA. It might be possible that the effective Hamiltonian in 



the extended phase space given by Eq. (43) cannot be applied to DW systems because 



it is derived with the assumed squeezed Gaussian wavepacket given by Eq. § in TDVA. 
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FIG. 9: (Color online) 3D plot of as functions of x and t of the symmetric DW model 

(xo = 0.0,^0 = 0.5). 

Nevertheless, quantum chaos in DW systems first pointed out Pattanayak and Schieve ^ by 
using the TDVA, is supported by our SA calculation, which is consistent with Ref. \W\. The 
GWM provides us with an efficient and physically-transparent calculation method, clarifying 
the relation between quantum and classical mechanics. The GWM is considered to be best 
applied to dynamics in HO and AO with a small anharmonicity. For a better description of 
quantum DW systems, it might be necessary to adopt extended GWMs with superimposed 
multiple Gaussian wavepackets (see Ref. [25], related references therein), which are much 
sophisticated and complicated than the original Heller's GWM [3j. 
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FIG. 10: (Color online) The x dependence of at various t of the SDW model calculated 

by the SA (bold solid curves) and the TDVA (solid curves) (xq = 0.0, — 0.5). 
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